
wo 00/32826 PCT/US99/28403 

-1- 

METHODS AND KITS FOR IDENTIFYING INDIVIDUALS AT RISK 
OF DEVELOPING OSTEOPOROSIS 

INTRODUCTION 

5 This invention was made in the course of research 

sponsored by the National Institutes of Health, The U,S, 
Government may have certain rights in this invention , 
BACKGROUTND OF THE INVENTION 

Osteoporosis is a disabling condition characterized by 

10 low bone density and fragile bones. The term 

"osteoporosis" , encompasses diseases of many different 
etiologies that all result in a reduction in the mass of 
bone per unit volume. Osteoporosis can be seen without 
association with other disease (idiopathic) and also in 

15 conjunction with other diseases (e.gr., rheumatoid 

arthritis, diabetes, calcium deficiency) . The decreased 
bone mass characteristic of the disease is not due to a 
decrease in the ratio of mineral to organic phase of bone, 
nor to any known abnormality in the mineral or organic 

20 matter of bone. Histological changes seen include a 

decrease in cortical thickness and a decrease in the number 
and size of the trabeculae of cancellous bone. 
Osteoporosis is the most common of all metabolic bone 
diseases and is seen commonly in the elderly. 

25 Since bone mass is decreased in osteoporosis, the rate 

of bone resorption exceeds that of bone formation. There 
is a normal age-related loss of bone mass in adults. 
However, this age-related loss of bone begins earlier and 
proceeds more rapidly in women. The reasons for the age- 

30 related bone loss are not known, although risk factors 
identified have included race (white women being at 
greatest risk) , sex (females) , sedentary lifestyle, and 
less muscular development. Other factors associated with 
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risk of developing osteoporosis include early menopause, 
smoking, and excessive alcohol consumption. 

Diagnosis of osteoporosis is most often done in 
conjunction with a study of bone density by radiography. 
5 Although clinical laboratory tests such as levels of 
calcium and phosphorus in blood can be examined, these 
measures are usually normal in osteoporotic patients. Only 
about 20% of postmenopausal women with osteoporosis exhibit 
hypercalciuria, or increased excretion of calcium in urine. 
0 Therefore, such laboratory findings are not indicative of 
the presence of disease, and clearly would not be 
indicative of risk of developing disease. To date, the 
prediction of risk of developing disease relies on family 
history of the disease. However, no genetic test is 
currently available to screen individuals. 

The ethnic differences shown to exist in the 
propensity to develop osteoporosis have led researchers to 
believe that genetic factors play a dominant role in the 
etiology of this disease. Several genes have been shown to 
be associated with low bone density and research has 
focused on identifying those genes that may act as markers 
of disease. Common allelic variations of the vitamin D 
receptor gene have been found to be associated with 
decreased bone density in certain populations, including 
premenopausal women and young girls (Wood, R.J. and Fleet, 
J.C. Ann, Rev. Nutrit. 1998 18:233-258). Bone mineral 
density has also been associated with genetic variation in 
the estrogen receptor gene, both by itself and in 
conjunction with variations in the vitamin D receptor gene 
(Willing et al . J. Bone Min. Res. 1998 13:695-705). In 
Japanese women, the HLA-A*24 -B*07-DRB*01 halotype has been 
linked to low peak bone mass (Tsuji et al . Hum. Irmunol . 
1998 59:243-249) . A variant of the gene encoding 
transforming growth factor-beta 1 has also been associated 
with low bone mass in osteoporotic women and with low bone 
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mass and increased bone turnover in normal women (Langdahl 
et al. Bone 1997 20:289-294), A polymorphism of the COLIAl 
gene has been identified as a potential marker for low bone 
mass and vertebral fracture in women (Grant et al . Nat. 
5 Genet. 1996 14:203-205). However, none of these studies 
has examined the relationship between tumor necrosis factor 
alpha receptor 2 gene and osteoporosis or bone mass in 
women . 



10 determined that there was a gene or genes on chromosome 1 
of humans that was linked to low bone density. However, 
there are many potential genes on this chromosome that 
could be involved in the etiology of low bone density. 



15 developing osteoporosis can be identified by assessing the 
genotype of the tumor necrosis factor alpha receptor gene 
in a sample of DNA obtained form the individual . A 
polymorphism in the gene on chromosome 1 for tumor necrosis 
factor alpha receptor 2 has now been shown to be associated 

2 0 with low bone density. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide a 
method of screening individuals for risk of developing 
osteoporosis. This method is based on assessment of the 
25 genotype of the tumor necrosis factor alpha gene in a DNA 
sample obtained from the individual . The presence of 
polymorphisms at nucleotides 593, 598 and 620 in exon 10 of 
the tumor necrosis factor alpha 2 receptor gene is 
associated with low bone density. 

3 0 Another object of the present invention is to provide 

a kit for screening individuals for risk of developing 
osteoporosis which comprises a means for assessing the 
genotype of the tumor necrosis factor alpha 2 receptor gene 
in an individual . 



Devoto et al. {Eur. J. Hum. Genet. 1998 6:151-157) 



It has now been found that individuals at risk of 
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DETAILED DESCRIPTION OF THE INVENTION 

In the present invention a method is provided for 
screening individuals for the presence of a polymorphism in 
the tumor necrosis factor alpha receptor 2 (TNFR2) gene. 
5 The presence of this genotype in individuals has now been 
shown to be associated with low bone density and the risk 
of developing osteoporosis. The association between TNFR2 
gene and low bone density was based on the observation that 
there is a polymorphism in this gene at a particular region 

10 of the gene sequence. This polymorphism was identified by 
Kaufman et al . (Hum. Genet. 1994 94:418-422), however, no 
characterization of its nature or examination of the 
relationship with bone density was disclosed. In the 
present invention, a relationship of this polymorphism with 

15 bone density has now been identified. 

In these experiments, two polymorphic sites within the 
TNFR-2 gene, including a microsatellite repeat within 
intron 4 and a three nucleotide variation in the 
3 ' untranslated region of the gene were assayed in three 

20 groups of individuals. Test group 1 was comprised of 

individuals from families used for linkage analysis of the 
low BMD (mean bone density) trait. Test group 2 was made 
up of unrelated individuals from an osteoporosis clinic. 
The third group was made up of unrelated individuals from a 

25 general clinic population with no testing for or knowledge 
of their bone status. This third group was used as 
reference for allele frequency comparisons. The 
composition and mean bone density Z-scores for test groups 
1 and 2 are shown in Table 1 , 
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Table 1 : Results of ANOVA among the males and females 

comprising the two test groups for age, body mass 
index and bone density 





GROUP 1 


GROUP 2 


P 


F ratio 




Male 








Value 


\ar ; 




Female 


Male 


Female 




Nuinber 


22 


44 


12 


81 






Age 


56.2 


57.3 


43.7 


52 . 6 


0 . Oil 


3 .846 


(yrs + SE) 


(2.5) 


(1.9) 


(3.8) 


(1.5) 




V J , ±z>^ } 


Body Mass 


28.4 


25 . 9 


25.0 


24 . 8 


0 . 084 


2 .259 


Index (+SE) 


(1.3) 


(0.9) 


(1.7) 


(0,6) 




(3, 149) 


Mean BHD Z- score {+SE) 




L2-L4 


-0.56 


-1.17 


-1.02 


-0.99 


0 .442 


0 . 901 




(0.38) 


(0.21) 


(0 .49) 


(0.14) 




(3,155) 


Femoral neck 


-0.50 


-0 . 74 


-0,73 


-0 .74 


0.790 


0.348 




(0 .27) 


(0.13) 


(0.39) 


(0.09) 




(3, 150) 


Trochanter 


0 .04 


-0.64 


-0 . 84 


-0.59 


0 . 054 


2 .611 




(0.31) 


(0.16) 


(0.30) 


(0.11) 




(3,146) 


Ward' s 


-0.21 


-0 . 76 


-0.61 


-0.70 


0.232 


1 .446 


triangle 


(0.22) 


(0.16) 


(0.45) 


(0.11) 




(3,146) 



15 Only age was significantly different among the groups 

(p=0.01l), due to the younger mean age of the males in the 
second group. However, ANOVA of age as the dependent 
variable with genotype for the two polymorphisms as class 
variable indicated that age was unaffected by genotype. 

2 0 Therefore the results of genotype testing in the two groups 
were pooled for subsequent analyses . 

The microsatellite repeat within intron 4 of the TNFR2 
gene has been reported previously to have 5 alleles 
(Beltinger et al . Genomics 1996 35:94-100). In these test 

25 groups, 6 alleles were observed with a statistically 

significant difference in the frequency of distribution 
between the two groups (p<0,025. Table 2) . A total of 18 
different genotypes would be expected from 6 alleles 
assuming Hardy -Weinberg equilibrium. The observed 

30 frequency distribution of genotypes of the intron 4 

microsatellite repeat for both the test and control groups 
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did not differ significantly from expected (Test group Chi- 
square=20 .54, df=19, N.S,; control group Chi-square=10 . 08 , 
df=19, N.S.). Comparison of the genotype distributions 
between the test and control groups was also not 
significant (Chi-square=22 , 46 , df^l9, N.S.). 

The polymorphism within the 3'UTR of the TNFR2 gene 
was previously characterized as having 4 alleles that were 
identified by different banding patterns under 
electrophoretic conditions that detect single strand 
conformation differences (Kaufman et al . Hum. Genet. 1994 
94:418-422; White et al , Genes, Chrow. , Cane. 1993 7:102- 
108; and Kaufman et al . Hum. Molec. Genet. 1993 2:824). 
Using a different conformation sensitive gel technique, 5 
alleles have now been resolved, each as a single band 
(Table 2) . The allele frequency distribution was compared 
between the test and control groups, with no statistically 
significant differences. The frequency distribution of the 
15 genotypes was in Hardy -Weinberg equilibrium in both . 
groups (Test group Chi-square= 8.44, df=14, N.S,; Control 
group Chi-square=22 . 57, df=14, N.S.). However, comparison 
of the genotype distributions for the 3*UTR polymorphism 
between the test and control groups was statistically 
significant (Chi -square=4 8 . 70 , df=14, p<0.001). 
Table 2: Allele frequency distribution for the two TNFR2 



polymorphisms in test group (n=159) compared with 
control group (n=141) 



Allele 


Intron 4 Polymorphism* 


UTR Polymorphism** 


Test group 


Controls 


Test group 


Controls 


1 


0.195 


0 . 156 


0 . 059 


0 . 038 


2 


0.037 


0 . 080 


0 .229 


0.240 


3 


0 . 573 


0.551 


0 .154 


0,137 
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4 


0 . 186 


0 . 193 


0 .374 


0 .427 


5 


0 . 006 


0 . 004 


0.182 


0.156 


6 


0.003 


0 . 015 


NA*** 


NA 



*Chi Square = 14.74, df=:5, p<0.025 



5 **Chi Square = 8.24, df=4, p<0,10 
***NA=Not applicable 

The nature of the UTR polymorphism was established by 
sequencing of PGR product from individuals who had been 
scored as homozygous for alleles 2, 3, 4 and 5 or 
10 heterozygous for combination of allele 1 and the other four 
alleles. No 1/1 homozygotes were observed. The 
polymorphism is actually a combination of three polymorphic 
sites within a span of 28 nucleotides, affecting 
nucleotides 593, 598 and 62 0 of exon 10 of the TNFR2 gene 
15 (Genbank Accession #U52165; SEQ ID N0:1) . The five 

observed alleles are defined as summarized in Table 3 . 
Table 3: Alleles of the TNFR2 UTR polymorphism as defined 
by nucleotide composition at three variable 
sites. Numbering is based on Genbank Accession 
20 #U52165 



Alleles 


Nucleotide 
593 


Nucleotide 
598 


Nucleotide 
620 


1 


A 


G 


T 


2 


A 


T 


T 


3 


G 


T 


C 


4 


G 


T 


T 


5 


A 


T 


C 



Analysis of variance (ANOVA) was performed on BMD 2- 
scores at the spine, femoral neck, trochanter and Ward's 
triangle with respect to genotype at each polymorphic site 
30 (Table 4) . 



wo 00/32826 PCT/US99/28403 

-8- 

Table 4: Results of ANOVA of BMD Z-scores with genotype 
for each of the two TNFR2 polymorphisms as the 
grouping variable in the study group (n=159) 





TNFR2 Polymorphism 




Intron 4 
Microsatellite 


3 'UTR 




P value 


F ratio 
(df) 


P value 


F ratio 
(df) 


L2-L4 


0 . 044 


1.937 
(10, 146) 


0 , 170 


1 . 386 
(13, 145) 


Femoral 
Neck 


0 . 120 


1.574 
(10, 141) 


0 . 140 


1 .459 
(13, 140) 


Trochanter 


0.199 


1.373 
(10, 137) 


0 . 127 


1 .494 
(13,136) 


Ward' s 
triangle 


0.199 


1.372 
(10, 137) 


0 . 149 


1.439 
(13,136) 



No statistically significant relationship was observed 
(i.e. p<0,025). Analysis of the genotype data using Tukey- 
Kramer post hoc pair-wise comparisons indicated that the 
mean Z- score for spinal BMD of individuals having genotype 
1/2 for the UTR polymorphism was significantly different 
from those with genotype 2/2. Analysis of genotype data 
for the CA polymorphism by the Tukey- Kramer test suggested 
no significant differences. However, when the UTR 
genotypes were placed in order with respect to mean spinal 
BMD Z-scores, it was found that the lowest mean spinal BMD 
occurred in the genotype groups that had allele 1, the 
rarest allele (Table 5) . 

Table 5: Genotypes for the TNFR2 UTR polymorphism with 



frequency and mean spinal or femoral neck BMD Z- 
scores for the study group 



GENOTYPE 


FREQUENCY 


SPINAL BMD 2- 
SCORE* 
(Mean + SE) 


FEMORAL NECK BMD 
Z- SCORE** 
(Mean + SE) 


1/2 


5 


-2 . 814 (0 .633) 


-0.727 (0.469) 
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.469) 
.663) 
.354) 
.215) 
.354) 
.184) 
.260) 
.419) 
.195) 
.354) 
.331) 
.221) 
.282) 



15 *ANOVA, p=0.17, R==0.11 
**ANOVA, p=0.14, R2=0.12 

The mean femoral neck BMD Z- scores are shown for 
comparison, and indicate a similar trend in that the two 
lowest mean Z-scores occur in groups having allele 1. 

20 Thus, these analysis show that the frequency of the 

rarest allele of this gene is associated with low bone 
density. Further, this association between the genotype of 
TNFR2 and bone density is statistically significant 
(p=0.008) . Accordingly, assessment of the genotype of the 

25 TNFR2 gene of an individual is useful in identifying 

individuals with low bone density at risk for osteoporosis. 

In the method of the present invention, a DNA sample 
is obtained from an individual. The genotype of the TNFR2 
gene is then assessed to determine if the individual is at 

30 risk for developing osteoporosis. For example, detection 
of a polymorphisms at nucleotides 593, 598 and 620 of exon 
10 of the TNFR2 gene (SEQ ID NO:l) is indicative of an 
individual at risk of developing osteoporosis. 
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1/5 



-1.718 (0.708) 



1.445 (0 



1/3 



-1.345 (1.00) 



1.620 (0 



1/4 



1.316 (0.500) 



-0.454 (0 



3/4 



21 



-1.225 (0.309) 



0.920 (0 



2/3 



2/4 



26 



•1.218 (0.535) 



-0.494 (0 



1.124 (0.277) 



-0.777 (0 



2/5 



13 



-0.980 (0.393) 



•0.342 (0 



3/3 



-0.778 (0.578) 



-0.458 (0 



4/4 



23 



0.776 (0.295) 



-0.868 (0 



10 



3/5 



-0.750 (0.535) 



-1.108 (0 



5/5 



-0.719 (0.500) 



-1.130 (0 



4/5 



18 



•0.684 (0.334) 



-0.508 (0 



2/2 
Total 



11 



159 



0.023 (0.427) 



•0.022 (0 
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The present invention also relates to diagnostic kits 
comprising a means for identifying individuals at risk of 
developing osteoporosis based on assessment of the genotype 
of the TNFR2 gene. In a preferred embodiment, such means 
5 are capable of detecting polymorphisms at nucleotides 593, 
598 and 620 of exon 10 of the TNFR2 gene (SEQ ID NO:l). 
Such means can be developed routinely by those skilled in 
the art based upon known techniques for genotype screening 
such as PGR and gel electrophoresis. In one embodiment, 

10 the kit comprises primers such as exemplified by SEQ ID NO: 
4 and SEQ ID NO: 5 described in Example 3 herein. However, 
as will be obvious to those of skill in the art upon this 
disclosure, additional means useful in assessing the TNFR2 
genotype of an individual can also be used. In a preferred 

15 embodiment, the kits of the present invention further 

comprise a blood collection device for obtaining a blood 
sample from an individual as the source of DNA for 
genotyping. Other sources of DNA obtained from a patient 
which can be genotyped via the kits and methods of the 

20 present invention include both tissue samples and other 
biological fluids . 

Using the method and/or kits of the present invention, 
one of skill can screen individuals known to have a family 
history of low bone density and/or individuals with 

25 identified risk factors for osteoporosis and low bone 
density including, but not limited to, post -menopausal 
women and the elderly. Via these methods and kits, 
individuals having polymorphisms in the TNFR2 gene, and in 
particular polymorphisms detected at nucleotide 593, 598 

3 0 and 62 0 of exon 10 of the TNFR2 gene (SEQ ID NO:l), can be 
identified as being at risk of developing osteoporosis. 

The following non- limiting examples are provided to 
further illustrate certain aspects of the instant 
invention. 
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EXAMPLES 



Example 1: 



Sxabjects 



10 



' 15 



20 



25 



The first group of test subjects was composed of 66 
individuals unrelated to each other. These individuals 
were either the founders of or had married into families 
with low BMD. The second test group of 93 individuals were 
also unrelated to each other and to the first group and had 
been referred to the McGill Bone Center by physicians for 
evaluation of bone status. The third group was a random 
general clinic population (n=141) , that served as a control 
for allele frequency distribution in these studies. 

Blood was collected from participants and DNA was 
isolated by standard methods. Bone densitometry was 
performed on individuals in the first two groups in 
accordance with procedures described by Devoto et al , Eur, 
J. Hum. Genet. 1998 6:151-157 and Spotila et al . Calcif. 
Tissue Int. 1996 59:235-237. 

Example 2: Microsatellite Repeat in Intron 4 

The PGR primers for detection of the (CA) 
microsatellite repeat described by Beltinger et al . 
{Genomics 1996 35:94-100) were TNFR2 -CAF : 5 ' - 

GTGATCTGCAAGATGAACTCAC-3 ■ (SEQ ID NO: 2) and TNFR2 -CAR : 5 ' - 
ACACCACGTCTGATGTTTCA-3 ' (SEQ ID NO: 3) , PGR conditions were 
denaturation at 94.5°C for 3 minutes; followed by 5 cycles 
of 94.5^*0 for 30 seconds, 65.0°C for 30 seconds, 72'*C for 30 
seconds, followed by 23 cycles of 94.5°C for 3 0 seconds, 
55.0°C for 3 0 seconds, 72.0°C for 3 0 seconds; followed by 
72.0''C for 5 minutes in a Perkin Elmer Thermocycler 9600. 
The primer TNFR2-CAF was labeled with y-^^P-ATP using T4 
polynucleotide kinase. Approximately 100-300 ng of DNA 
template was used with reaction components under conditions 
recommended by the manufacturer (Applied Biosystems) . 
Example 3: Polymorphism of 3» Untranslated Region 

PGR primers for detection of the polymorphism within 



# 



• 
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10 



15 



20 



25 



30 



the 3'UTR of TNFR2 as described by Kaufman et al , (Hum. 
Genet. 1994 94:418-422) were TNFR2-7:5'- 
AGGACTCTGAGGCTCTTTCT- 3 ' ( SEQ ID NO : 4 ) and TNFR2 -8:5'- 
TCACAGAGAGTCAGGGACTT-3 ' (SEQ ID NO: 5). PGR conditions 
were denaturation at 94.5*^0 for 3 minutes; followed by 28 
cycles of 94.5**C for 30 seconds, 62.0°C for 30 seconds, 
72.0°C for 3 0 seconds; followed by 72.0''C for 5 minutes. 
Reagents and instrumentation were the same as described for 
Example 2 , 

Using conformation sensitive gel electrophoresis 
(CSGE) , 5 differentially migrating fragments of this PGR 
product were distinguished. Gel conditions and product 
suppliers were similar to those described by Ganguly et al. 
{Proc. Natl Acad, Sci. USA 1996 93:13721-13725) with the 
following modifications. Tris- taurine -EDTA running buffer 
was made as a 20X stock with 1.78 M Tris, 0.58 M Taurine, 
4.0 mM EDTA, pH 9 . 0 and used at 0.5X. Formamide was 
deionized in batches and frozen in aliquots sufficient for 
a single gel. Sample loading buffer was prepared each time 
with 4 parts deionized formamide and 1 part loading buffer 
(5X=0.1% xylene cyanol, 0.1% bromophenol blue, 10 mM EDTA, 
95% formamide) . Standard sequencing gel plates and 
apparatus were used with 0,4 mm spacers . Samples were 
mixed with an equal volume of sample loading buffer, 
denatured for 5 minutes at 95°C, transferred to ice and 1,5 
^1 was immediately loaded. Gels were run at 350 V at room 
temperature for 16 to 18 hours, then dried and exposed to 
film. Resolution of 5 alleles was possible only when these 
conditions were adhered to. 
Example 4: Sequence Analysis 

Selected PGR fragments were sequenced according to 
manufacturer's recommendations using ABI 377 
instrumentation and Big-Dye chemistry (Applied Biosystems, 
Inc. Foster City, CA) in order to determine the basis for 
the 5 alleles. 



# 
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Exaznple 5: 



Statistical Analysis 



The two test groups were compared by ANOVA for age, 
body mass index and BMD Z-score at the spine, femoral neck, 
trochanter and Ward's triangle, and then combined for 
5 additional analyses. The frequency distribution of alleles 
for the two polymorphisms within the TNFR2 gene was 
analyzed by Chi Square. Association of bone density at the 
spine (L2-L4), femoral neck, trochanter and Ward's triangle 
with genotype for each of the two polymorphisms was 

10 evaluated by ANOVA. The BMD Z-scores at each of the 
anatomical sites were used as the independent variable 
since they were already corrected for age, gender and 
weight. Statistical significance was achieved at p<0.025, 
because of multiple testing. Pairwise comparisons of mean 

15 BMD Z-scores for the different genotypes were made for all 
pairs using the Tukey-Kramer post hoc test. Power analysis 
was performed to determine the lowest sample number 
required to observe a significant (p=0.025) effect of 
genotype on BMD. The statistical software package employed 

20 was JMP (Version 3.1, SAS Institute, Inc. Gary, N.C.), 



